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ABSTRACT 

In order to understand certain observed features of arc-like giant radio relics 
such as the rareness, uniform surface brightness, and curved integrated spec¬ 
tra, we explore a diffusive shock acceleration (DSA) model for radio relics in 
which a spherical shock impinges on a magnetized cloud containing fossil rela¬ 
tivistic electrons. Toward this end, we perform DSA simulations of spherical 
shocks with the parameters relevant for the Sausage radio relic in cluster CIZA 
J2242.8-1-5301, and calculate the ensuing radio synchrotron emission from re¬ 
accelerated electrons. Three types of fossil electron populations are considered: 
a delta-function like population with the shock injection momentum, a power-law 
distribution, and a power-law with an exponential cutoff. The surface brightness 
prohle of radio-emitting postshock region and the volume-integrated radio spec¬ 
trum are calculated and compared with observations. We hnd that the observed 
width of the Sausage relic can be explained reasonably well by shocks with speed 
Us ~ 3 X 10^ km s“^ and sonic Mach number Mg ~ 3. These shocks produce 
curved radio spectra that steepen gradually over (0.1 — lOji^br with break fre¬ 
quency z/br ~ 1 GHz, if the duration of electron acceleration is ~ 60 — 80 Myr. 
However, the abrupt increase of spectral index above ~ 1.5 GHz observed in 
the Sausage relic seems to indicate that additional physical processes, other than 
radiative losses, operate for electrons with 7 e > 10^. 

Subject headings: acceleration of particles — cosmic rays — galaxies: clusters: 
general — shock waves 
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1. INTRODUCTION 


Radio relics are diffuse radio sources found in the outskirts of galaxy clusters and they 
are thought to trace synchrotron-emitting cosmic-ra y (CR) electrons accelerated via diffu¬ 


sive shock acceleration ( pSA) at cluster shoc ks fe.g. lEnsslin et alJll998l: iBagchi et al.ll2006 
van Weeren et al.ll2010l: IBriiggen et al.N2012l) . So far several dozens of clusters have been 


observed to have radio relics with a variety of morphologies and most of them are consid¬ 


ered to be associated w ith cluster merger activities (see for reviews, e.g., iFeretti et al.l 12012 


Brunetti fc Jonesiboid ). For instance, double radio relics, such as the ones i n ZwC10008.8-|-521 5, 


are th o ught to reveal the bow s hocks induced by a binary major merger (Ivan Weeren et ah 


201 Ibl: Ide Gasperin et al.l 120141) . On the other hand, recently it was shown that shocks in¬ 


duced by the infall of the warm-hot intergalactic medium (WHIM) along adjacent hlaments 
into the hot intracluster medium (ICM) can efficiently accelerate CR elec trons, and so the y 


could be responsible for some radio relics in the cluster outskirts fse e. e.g..lHong et al 


The radio r elic 1253-1-275 in Coma c luster observed in both radio flBrown fc Rudni eld 120111) 


2014Jj 


and X-ray fjOgrean fc Bruggenll2013l) provides an example of such infall shocks. 


The so-called Sausage relic in CIZA J2242. 8-1-5301 {z = 0.192) contains a thin arc-like 
structure of ~ 55 kpc width and ~ 2 Mpe length, which could be represented by a portion 
of spherical shell with radius ~ 1.5 Mpe ( van Weeren et al. 201ol) . Unique features of this 
giant radio relic include the nearly uniform surface brightness along the length of the relic 
and t he strong polarization o f up to 50 — 60% with magnetic held vectors aligned with the 
relic (Ivan Weeren et al.l 120101) . A temperature jump across the reli c that corresponds to a 


M, 2.54 — 3.15 shock has been detected in X-ray observations flAkamatsu fc Kawahara 


20131: lOgrean et al.l 120141) . This was smaller than Mg ~ 4.6 estimated from the above radio 


observation. Several examples of Mpc-scale rad io relics include the To othbrush relic in 
IRXS J 0603.3 with a peculiar linear morph ology (Ivan Weeren et al.ll2012h and the relics in 
A3667 f Rottgering et al. 1997 ) and A3376 ( Bagchi et al. 2006 ). The shock Mach numbers 
of radio relics estimated based on X-ray observation are often lower than those inferred 
from the radio spec tral index using the DSA model, for instance, in the Toothbrush relic 
fjOgrean et al.l 120131) and in the radio relic in A2256 flTrasatti et al.l 120151) . Although such 
giant radio relics are quite rare, the fr action of X-ray lum inous clusters hosting some radio 
relics is estimated to be ~ 10 % or so (IFeretti et al.ll2012l) . 


Through a number of studies using cosmological hydrodynamical simulations, it has 
been demonstrated that during the process of hierarchical structure formation, abundant 
shocks are produced in the large-scale structure of the universe, especially in clusters (e.g.. 


Rvu et al. 

20031: 

Pfrommer et al. 2 

i006; 

Skillman et al. 

2008; 

Hoeft et al. 

2008: 

Vazza et al. 

2009; 

Vazza et al 

2011; 

Hong et al 

12014). Considering that the characteristic time-scale of 
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cluster dynamics including mergers is tdyn ~ 1 Gyr, typical cluster shocks are expected to 
last for about the same period. Yet, the number of observed radio relics, which is thought 
to trace such shocks, is still limited. So it is plausible to conjecture that cluster shocks may 
‘turn on’ to emit synchrotron radiation only for a fraction of their lifetime. One feasible 
scenario is that a cluster shock lights up in radio when it sweeps up a fossil cloud, i.e., a 
magnetized ICM gas with fossil relativistic electrons left over from either a radio jet from 
AGN or a previous episode of shock/turbulence acceleration (see the discussions in Section 
2.5 and Figure 1). 


Pre-exiting seed electrons and/or enhanced magnetic helds are the requisite conditions 
for possible lighting-up of relic shocks. In particular, the elongated shape with uniform 
surface brightness and high polarization fraction of radio emission in the Sausage relic, may 
be explained, if a Mpc-scale thermal gas cloud, containing fossil relativistic electrons and 
permeated with regular magnetic field of a few to several fiG, is adopted. A more detailed 
description will be given later in Section 2.5. In this picture, fossil electrons are expected 
to be re-accelerated for less than cloud-crossing time (< Rdoud/us ~ 100 Myr), which is 
much shorter than the cluster dynamical time-scale. In addition, only occasional encounters 
with fossil clouds combined with the short acceleration duration could alleviate the strong 
constraints on the PSA theory based on non-detectio n of 7 -ray emission from clusters by 
Fe rmi-LAT flAckerrnann et al.ll2014j: IVazza et al.ll2015[) . A similar idea has been brought up 
bv IShimwell et al.l (120151), who reported the discovery of a Mpc-scale, elongated relic in the 
Bullet cluster IE 0657-55.8. They also proposed that the arc-like shape of uniform surface 
brightness in some radio relics may trace the underlying reg ions of pre-existing, seed el e ctrons 
remaining from old radio lobes. On the other hand, lEnsslin fc Gopal-Krishnal fj200l[) 
suggested that radio relics could be explained by revival of fossil ra dio plasma by comp r ession 
due to a passage of a shock, rather than DSA. In a follow-up study, lEnsslin fc Briiggenl (1200211 
showed using MHD simulations that a cocoon of hot radio plasma swept by a shock turns into 
a hlamentary or toroidal structure. Although this scenario remains to be a viable explanation 
for some radio relics, it may not account for the uniform arc-like morphology of the Sausage 
relic. 


It is now well established, through observations of radio halos/relics and Faraday ro¬ 
tation measures of background radio sources, that the IGM is permeated with /iG-level 
magnetic fields (e.g. iBonafede et al.l 1201 ll: iFeretti et al.l 1201211 . The observed radial profile 
of magnetic field strength tends to p eak at the center with a few /iG and decrease outward 
to ~ 0.1/iG in the cluster outskirts (IBonafede et al.l 1201011 . A variety of physical processes 
that could generate and amplify magnetic fields in the IGM have been suggested: primor¬ 
dial processes, plasma processes at the recombination epoch, and Biermann battery mecha¬ 
nism, combined with turbulence dynamo, in addition to galactic winds and AGN jets (e.g. 
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Rvu et al.l 120081: IPolag et al.l I2OO8I: iBriiggen et al.l I2OI2I: IRvu et al.l I2OI2I: I Choi l2014l) . Given 
the fact that ~ 5 /iG helds are required to e xplain the amplitude and width of the obs erved 


radio flux prohle of the Sausage relic (e.g. Ivan Weereu et al.l I 2 OIOI: iKaug et al.l l2012l) . the 


presence of a cloud with enhanced magnetic helds of several fiG mig ht be preferred to the 
backg round helds of ~ O.l/rG in the cluster periphery. Alternatively, llapichino Sz Briiggen 
fl2012l) showed that the postshock magnetic helds can be amplihed to ~ 5 — 7/iG level lead¬ 
ing to high degrees of polarization, if there exists dynamically signihcant turbulence in the 
upstream region of a curved shock. Although it is well accepted that magnetic helds can 
be amplihed via various plasma instabilities at collisionless shocks, the dependence on the 
shock parameters such as the shock sonic and Alfvenic Mach numb ers, and the obliqu ity of 
background magnetic helds remains to be further investigated fsee lSchure et al.ll2012l) . For 
example, the acceleration of protons and ensuing magnetic held amplihcation via resonant 
and non-resonan t streaming instabilities are f ound to be inehective at perpendicular shocks 
( Bell 1978 . 2004 : Gaprioli fc Spitkovskv 2014). 


In several studies using cosmological hydrodynamical simulations, synthetic radio maps 
of simulated clusters were constructed by identifying shocks and adopting models for PSA of 


electrons and magnetic held amplihcation fiNuza et ah 
2 OI 3 I: iHong et al.ll2015l) . In particular, IVazza et ah 


2OI2I: 


Vazza et ^l2012a : Skillman et al 
( 2012bl) demonstrated, by generating 


mock radio maps of simulated cluster samples, that radio emission tends to increase toward 
the cluster periphery and peak around 0.2 — O.SRvir (where Rvir is the virial radius), mainly 
because the kinetic energy dissipated at shocks peaks around 0.2i?vir- As a result, radio 
relics are rarely found in the cluster central regions. Re-acceleration of fossil relativistic elec- 
trons by cosmologic al shocks during the large scale structure formation has been explored by 
Pinzke et al.l (120131 1. The radio emitting shocks in these studies look like segments of spheri¬ 
cal shocks, moving from the cluster core region into the periphery. We presume that they are 
generated mostly as a consequence of major mergers or energetic infalls of the WHIM along 
adjacent hlaments. So it seems necessary to study spherical shocks propagating through the 
cluster periphery, rather than interpreting the radio spec tra by PSA a t steady planar shocks, 
in order to better understand the nature of radio relics (jKane 2015a ,b). 


According to the PSA theory, in the case of a steady planar shock with constant post¬ 
shock magnetic held, the electron distribution function at the shock location becomes a 
power-law of fe{p, Ts) oc p~'^, and so the synchrotron emissivity from those electrons becomes 
a power-law of ju{rs) oc The power-low slopes depend only on the shock sonic Mach 

number. Mg, and are given as q = — 1) a nd + 3)/2(M^ — 1) for the 


gasdynamic shock w ith the adiabatic index 7 ^ = 5/3 (jPrurvlll983l : lBlandford fc Eichlerll 19871: 
Ensslin et al.lll998l ). Here we refer ctinj as the injection spectral index for a steady planar 
shock with constant postshock magnetic held. Then, the volume-integrated synchrotron 
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spectrum downstream of the shock also becomes a simple power-law oi Jy = f jy(r)dV oc 
with the spectral index Ay = Oinj + 0.5 above the break frequency, z/br, s ince electrons 


cool via synchrot ron and inverse-Compton (IC) losses behind the shock (e.g.. lEnsslin et ah 


1998 


Kang 2011 )^1 Such predictions of the DSA theory have been applied to explain the 


observed properties of radio relics, e.g., the relation between the injection spectral index 
and the volume-integrated spectral index, and the gradual steepening of spatially resolved 
spectrum downstream of the shock. 


Kang et al.l (120121 1 performed time-dependent, DSA simulations of CR electrons for 


steady planar shocks with Mg = 2 — 4.5 and constant postshock magnetic helds. Several mod¬ 
els with thermal leakage injection or pre-existing electrons were considered in order to repro¬ 
duce the surface brightness and spectral aging prohles of radio relics in CIZA J2242.8-1-5301 
and ZwC10008.8-1-5215. Adopting the same geometrical structure of radio-emitting volume 
as described in Section 2.5, they showed that the synchrotron emission from shock acceler¬ 
ated electrons could explain the observed prohles of the radio hux, Sy{R), of the Sausage 
relic, and the observed prohles of both Sy{R) and ay{R) of the relic in ZwC10008.8-1-5215. 
Here R is the distance behind the projected shock edge in the plane of the sky. 

In the case of spherically expanding shocks with varying speeds and/or nonuniform 
magnetic held prohles, on the other hand, the electron spectrum and t he ensuing ra dio 
spectrum could deviate from those simple power-law forms, as shown in iKangl fl2015al]bl) . 
Then even the injection slope should vary with the frequency, i.e., ainj(z^). Here we follow the 
evolution of a spherical shock expanding outward in the cluster outskirts with a decreasing 
density prohle, which may lead to a curvature in both the injected spectrum and the volume- 
integrated spectrum. Moreover, if the shock is relatively young or the electron acceleration 
duration is short (< 100 Myr), then the break frequency falls in z/br ~ 1 GHz and the 
volume-integrated spectrum of a radio relic w ould steepen gradually with the spectral index 


from ttinj to Oinj + 0.5 over (0.1 — 10)z/br (e.g. lKangil2015bl) . 


In the case of the Sausage relic, Ivan Weeren et al.l (120101) and IStroe et ahl (120131) orig¬ 


inally reported observations of a 


jnj 


0.6 and a = Ay 


1.06, which imply a shock 


of Mg ^ 4.6. IStroe et ahl (j2014bl) . however, found a spectral steepening of the volume- 


integrated spectrum at 1 6 GHz, whic h would be inconsistent with the DSA model for a steady 
planar shock. Moreover, IStroe et ahl (l2014al) . by performing a spatially-resolved spectral £t- 
ting, revised the injection index to a steeper value, Oinj ~ 0.77. Then, the corresponding 


^Note that radio observers commonly use ‘a’ as the spectral index of the flux density, Sy oc z/““ for 
unresolved sources, so in that case a is the same as Ay. Here ay{r) is defined as the spectral index of the 
local emissivity, jy (r). See Equations dUl-dini). 
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shock Mach number is reduced to Mg ^ 2.9. They also suggested that the spectral age, cal¬ 
culated under the assumption of freely-aging electrons downstream of a steady planar shock, 
mig ht not be compatible w ith the shock speed estimated from X-ray and radio observations. 
Also iTrasatti et al.l (120151) reported that for the relic in A2256, the volume-integrated index 
steepens from ^ 0.85 for z/ = 351 — 1369 MHz to Ai, ^ 1.0 for z/ = 1.37 — 10.45 GHz, 
which was interpreted as a broken power-law. 


Discoveries of radio relic shocks with Mg 


3 in recent years have brought up the 


need for more accurate understanding of injections of protons and electrons at weak colli 


sionle ss shocks, especially at high plasma beta {(3p ~ 50 — 100) ICM plasmas (e.g. iKang et ah 
20141) . Here /dp is the ratio of the gas to magnetic held pressure. Injection of electrons 


into the Fermi Ist-order process has been one of long-standing problems in the DSA theory 
for astrophysical shocks, because it involves complex plasma kine tic processes that can be 


studied only through full Par ticle-in-Cell (PIC) simulations (e.g. lAmano fc Hoshinol 12009 


Riquelme fc Spitkovskvi 1201111 . It is thought that electrons must be pre-accelerated from 


their thermal momentum to several times the postshock thermal proton momentum to take 
part in the DSA process, and electron injection is much less efficient than proton injec¬ 
tion due to smaller rigidity of electrons. Several recent studies using PIC simulations have 
shown that some of incoming protons and electrons gain energies via shock drift acceleration 
(SDA) while drifting along the shock surface, and then the particles are reflected toward the 
upstream region. Those reflected particles can be scattered back to the shock by plasma 


power-law suprathermal populations (e.g.. 


Guo et ah 

2014a 

b; 

Park et ah 

2015) 


Such ‘self 


pre-acceleration’ of thermal electrons in the foreshock region could be sufficient enough even 
at weak shocks in high beta ICM plasmas to explain the observed flux level of radio relics. 
In these PIC simulations, however, subsequent acceleration of suprathermal electrons into 
full DSA regime has not been explored yet, because extreme computational resources are 
required to follow the simulations for a large dynamic range of particle energy. 


The main reasons that we implement the fossil electron distribution, instead of the shock 
injection only case, are (1) the relative scarcity of radio relics compared to the abundance of 
shocks expected to form in the ICM, (2) the peculiar uniformity of the surface brightness of 
the Sausage relic, and (3) curved integrated spectra often found in some radio relics, implying 
the acceleration duration < 100 Myr, much shorter than the cluster dynamical time. 

In this paper, we consider a DSA model for radio relics; a spherical shock moves into 
a magnetized gas cloud containing fossil relativistic electrons, while propagating through 
a density gradient in the cluster outskirts. Specihcally, we perform time-dependent DSA 
simulations for several spherical shock models with the parameters relevant for the Sausage 






















7 


relic. We then calculate the surface brightness prohle, 1,^, and the volume-integrated radio 
spectrum, J^, by adopting a specihc geometrical structure of shock surface, and compare 
them with the observational data of the Sausage relic. 

In Section 2, the DSA simulations and the model parameters are described. The com¬ 
parison of our results with observations is discussed in Section 3. A brief summary is given 
in Section 4. 


2. DSA SIMULATIONS OF CR ELECTRONS 


2.1. ID Spherical CRASH Code 

We follow the evolution of the CR electron population by solving the following diffusion- 
convection equation in the one-dimensional (ID) spherical geometry: 


dge dge _ 

dt dr 3r^ dr 


1 d{r‘^u) f dg, 


dy 


~ %£ + 


2 


r‘^D{r,p) 


dge 

dr 


p 


A 

dy \p 




( 1 ) 


where ge{r,p,t) = fe{r,p,t)p^ is the pitch-angle-averaged phase space distribution function 
of electrons, u{r,t) i s the flow vel ocity and y = ln(p/mec) with the electron mass mg and 
the speed of light c fjSkillind 119751) . The spatial diffusion coefficient, D{r,p), is assumed to 
have a Bohm-like dependence on the rigidity. 


,2„-l 


D{r,p) = 1.7 X lO^^cm^s 


R(r) 
1 /iG 


-1 


P 

rripC 


( 2 ) 


The cooling coefficient h{p) = —dp/dt accounts for radiative cooling, and the cooling time 
scale is dehned as 


p 


Rad(7e) = -rr\ = 9-8 X 10 
h[p) 


Be 


5 pG 


-2 


f—r' 

VlOV 


(3) 


where 7 e is the Lorentz factor of electrons. Here the ‘effective’ magnetic held strength, 
+ B^^^ with Rrad = 3.24 /iG(l -|- takes account for the IC loss due to the 
cosmic background radiation as well as synchrotron loss. The redshift of the cluster CIZA 
J2242.8+5301 is z = 0.192. 

Assuming that the test-particle limit is applied at weak cluster shocks with Mg < 
several (see Table 1), the usual gasdynamic conservation equations are solved to follow the 
background how s peed, nfr, t), using th e ID spherical version of the CRASH (Cosmic-Ray 
Amr SHock) code (IKang fc JonesI 120061) . The structure and evolution of u{r,t) are fed into 















Equation ([T]), while the gas pressure, Pg{r,t), is used in modeling the postshock magnetic 
held prohle (see Section 2.3). We do not consider the acceleration of CR protons in this study, 
since the synchrotron emission from CR electrons is our main interest and the dynamical 
feedback from t he CR proton press ure can be ignored at weak shocks {Mg < 3) in the test- 
particle regime (iKang fc Rvull20131 1. In order to optimize the shock tracking of the CRASH 
code, a comoving frame that expands with t he instantaneo us shock speed is adopted. Further 
details of DSA simulations can be found in iKand (j2015al) . 


2.2. Shock Parameters 


To set up the initial shock structure, we adopt a Sedov self-similar blast wave prop¬ 
agating into a uniform static medium, which can be specihed by two parameters, typi¬ 


cally, the explosion en ergy, Eq, and the background density, po flOstriker fc McKeel Il988 


Rvu fc Vishniad Il99l[) . For our DSA simulations, we choose the initial shock radius and 
speed, j and Ug^i, respectively, and adopt the self-similar profiles of the gas density p(r), 
the gas pressure Pg{r), and the flow speed u{r) behind the shock in the upstream rest-frame. 
For the hducial case (SAl model in Table 1), for example, the initial shock parameters are 
Tg^i = 1.3 Mpc, and Ug^i = 3.3 x 10^ km s“^. For the model parameters for the shock and 
upstream conditions, refer to Table 1 and Section 3.1. 

We suppose that at the onset of the DSA simulations, this initial shock propagates 
through the ICM with the gas density gradient described by a power law of r. Typical 
X-ray brightness profiles of observed clusters can be represented approx imately by th e so-call 


beta model for isothermal ICMs, p{r) oc [l-|-(r/rc)^] with ft ^ 2/3 f Sarazinlll986[) . In the 


outskirts of clusters, well outside of the core radius (r 3> rf), it asymptotes as p(r) oc r~^^. 
We take the upstream gas density of pup oc as the hducial case (SAl), but also consider 
Pup oc r~‘^ (SA3) and pup = constant (SA4) for comparison. The shock speed and Mach 
number decrea s e in tim e as the spherical shock expands, depending on the upstream density 
prohle. IKand f 2015bl) demonstrated that the shock decelerates approximately as Ug oc t~^^^ 


for Pup = constant and as Ug oc t~^^^ for pup oc r~^, while the shock speed is almost constant 
in the case of pup oc r~^. As a result, nonlinear deviations from the DSA predictions for 
steady planar shocks are expected to become the strongest in SA4 model, while the weakest 
in SA3 model. In the hducial SAl model, which is the most realistic among the three 
models, the ehects of the evolving spherical shock are expected to be moderate. 


The ICM temperature is set as kTi = 3.35keV, adopted from lOgrean et all (120141) . 
in most of the models. Hereafter, the subscripts “1” and “2” are used to indicate the 
quantities immediately upstream and downstream of the shock, respectively. Although the 























9 


ICM temperature is known to decrease slightly in the cluster outskirt, it is assumed to 
be isothermal, since the shock typically travels only 0.2 — 0.3 Mpc for the duration of our 
simulations < 100 Myr. 


2.3. Models for Magnetic Fields 


Magnetic helds in the downstream region of the shock are the key ingredient that governs 
the synchrotron cooling and emission of CR electrons in our models. We assume that the 
fossil cloud is magnetized to /xG level. As discussed in the Introduction, observations indicate 
that the magnetic held strength decreases from ^ 1 ~ 10 in the core reg ion to ~ 0.1 — 1 /xG 
in the periphery of clusters fe.g.. iBonafede et al.ll2010l: iFeretti et al.ll2012ll. This co rresponds 
to the plasma beta of /dp ~ 50 — 100 in typical IGMs fe.g.. IRvu et al.l 2008 . 2012). On the 
other hand, it is well established that magnetic helds can be amplihed via resonant and 
non-resonan t instabilities induced by GR protons streaming upstream of strong shocks (IBell 


19781. 1200411. In addition, magne tic helds can be amplihed by turbulent motions behind 


shocks flGiacalone &: Jokipiil 120071) . Recent hybrid plasma simulations have shown that the 
magnetic held amplihcation factor due to streaming GR protons scales with the Alfvenic 
Mach number, Ma, and the C R proton acceleration efficiency as {6B/B)^ 3Ma{Pct, 2/Piu'i) 
flCaprioli fc Spitkovskvl 120141) . Here, 6B is the turbulent magnetic held perpendicular to the 
mean background magnetic held, Pcr ,2 is the downstream CR pressure, and is the 
upstream ram pressure. For typical radio relic shocks, the sonic and Alfv enic Mach number s 
are expected to range 2 < Mg < 5 and 10 < Ma < 25, respectively (e.g., IHong et al.ll2014l) . 


The magnetic held amplihcation in both the upstream and downstream of weak shocks is 
not yet fully understood, especially in high beta ICM plasmas. So we consider simple models 
for the postshock magnetic helds. For the hducial case, we assume that the magnetic held 
strength across the shock transition is increased by compression of the two perpendicular 
components: 

B2 (t) = B, v/1/3 + 2a{ty/3, (4) 

where Bi and B 2 are the magnetic held strengths immediately upstream and downstream 
of the shock, respectively, and a{t) = P 2 /P 1 is the time-varying compression ratio across the 
shock. For the downstream region (r < r^), the magnetic held strength is assumed to scale 
with the gas pressure: 

Bdn(r,() = B^(t} . |P,(r,^)/P,J(^)|■''^ (5) 

where Pg, 2 {t) is the gas pressure immediately behind the shock. This assumes that the ratio 
of the magnetic to thermal energy density is constant downstream of the shock. Since Pg{r) 
decreases behind the spherical blast wave, B^nir) also decreases downstream as illustrated 
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in Figure 2. This fiducial magnetic field model is adopted in most of the models described 
in Table 1, except SA2 and SA4 models. The range of B 2 {t) is shown for the acceleration 
duration of 0 < tage < 60 Myr in Table 1, reflecting the decrease of shock compression ratio 
during the period, but the change is small. 


In the second, more simplified (but somewhat unrealistic) model, it is assumed that 
Bi = 2 yuG and B 2 = 7 /iG, and the downstream magnetic fie l d stre ngth is constant, i.e., , 
B^n = B 2 for r < Tg. This model was adopted in iKang et ahl (120121) and also for SA2 and 
SA4 models for comparison in this study. 


Kangl fj2015bl) showed that the postshock synchrotron emission increases downstream 


away from the shock in the case of a decelerating shock, because the shock is stronger at 
earlier time. But such nonlinear signatures become less distinct in the model with decreasing 
downstream magnetic fields, compared with the model with constant downstream magnetic 
fields, because the contribution of synchrotron emission from further downstream region 
becomes weaker. 


2.4. Injection Momentum 


As described in the Introduction, the new picture of particle injection emerged from 
recent PIG simulations is quite different from the ‘clas sical’ thermal leak age injection model 
commonly employed in previous studies of DSA (e.g., iKang et ahl 120021) . However, the re¬ 
quirement of p > 3pth,p for particles to take part in the full Fermi Ist-order process, scattering 
back and forth diffusive ly across the slmck transition zone with thickness, A/ghock rgiPth,p), 
seems to remain valid fjCaprioli et al.ll2015l: IPark et ahl 120151) . Here, pth,p = \/‘2mpkBT2 is 


the most probable momentum of thermal protons with postshock temperature T 2 and Vg is 
the gyroradius of particles. In other words, only suprathermal particles with the gyro-radius 
greater than the shock thickness are expected to cross the shock transition layer. Hence, we 
adopt the traditional phenomenological model in which only particles above the injection 
momentum, pinj ~ 5.3 mpU,J a, are allowed to g et injected into the GR populations at the 
lowest momentum boundary (IKang et al.N2002l l. This can be translated to the electron 
Lorentz factor, 7e,inj = Pinj/'^eC ~ 30(ms/3000 km s“^)(3.0/ct). In the case of expanding 
shocks considered in this study, Pinj(t) decreases as the shock slows down in time. 
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2.5. Fossil Electrons 


As mentioned in the Introduction, one peculiar feature of the Sausage relic is the uni¬ 
form surface brightness along the Mpc-scale arc -like shape, which require s a special geomet¬ 
rical distribution of shock-accelerated electrons (Ivan Weeren et alJl 2011 al) . Some of previous 
studies adopted the ribbon-like curved shock surface and the downs tream swept-up volume 


viewed edge-on with the view ing extension angle if ~ 10 ° (e.g., Ivan Weeren et ahM 2010 
Kang et al.l 120121: iKangl l2015bf) . We suggest a picture where a spherical shock of radius 
Vs ~ 1.5 Mpc passes through an elongated cloud with width tCcioud ~ 260 kpc and length 
^cioud ~ 2 Mpc, hlled with fossil electrons. Then the shock surface penetrated into the cloud 
becomes a ribbon-like patch, distributed on a sphere with radius ~ 1.5 Mpc with the 
angle = 360° • Wcioud/ i^nrs) 10°. The downstream volume of radio-emitting, reaccel- 
er ated e l ectrons has the width A/( 7 e) ~ {ug/cr) ■ min[fage, ^rad( 7 e)], as shown in Figure 1 


of Kang f 2015bh . Hereafter the ‘acceleration age’, fage, is dehned as the duration of elec¬ 


tron acceleration since the shock encounters the cloud. This model is expected to produce 
a uniform surface brightness along the relic length. Moreover, if the acceleration age is 
^age ^ 100 M yr, the volum e-integrated radio spectrum is expected to steepen gradually over 


0.1-10 GHz flKang|l2015bll . 


There are several possible origins for such clouds of relativistic electrons in the ICMs: 
(1) old remnants of radio jets from AGNs, ( 2 ) electron populations that were accelerated 
by previous shocks and have cooled down below 7 e < 10^, and (3) electron populations 
that were accelerated by turbulence during merger activities. Although an AGN jet would 
have a hollow, cocoon-like shape initially, it may turn into a hlled, cylindrical shape through 
diffusion, turbulent mixing, or contraction of relativistic plasmas, as the electrons cool ra- 
diatively. During such evolution relativistic electrons could be mixed with the IGM gas 
within the cloud. We assume that the cloud is composed of the thermal gas of 7 ^ = 5/3, 
whose properties (density and temperature) are the same as the surrounding IGM gas, and 
an additional population of fossil electrons with dynamically insignihcant pressure. In that 
regard, our fossil electron cloud is different from hot bubbles considered in previous stu dies 
of the interaction of a shock with a hot rarehed bubble (e.g., lEns.s]in fc Bruggenll2002l) . In 
the other two cases where e lectrons were accelerated either by shocks or by turbulence (see. 


Farnsworth et al.l 120131) . it is natural to assume that the cloud medium should contain 


e.g. 

both thermal gas and fossil electrons. 


Three different spectra for fossil electron populations are considered. In the hrst hducial 
case {e.g., SAl model), nonthermal electrons have the momentum around pinj, which corre¬ 
sponds to 7 e,inj ~ 20 — 30 for the model shock parameters considered here. So in this model, 
seed electrons with ~ ye^nj are injected from the fossil population and re-accelerated into 
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radio-emitting CR electrons with 7 e > 10^. Since we compare the surface brightness prohles 
in arbitrary units here, we do not concern about the normalization of the fossil population. 


In the second model, the fossil electrons are assumed to have a power-law distribution 
extending up to 7e ;§> 10^, 

/e,up(p) = /o ■ f—) ■ (6) 

\Pinj/ 

For the modeling of the Sausage relic, the value of s is chosen as s = 2ainj -l- 3 = 4.2 with 
ctinj = 0.6 (SAlp model). 


As mentioned in the Introduction, the volume-integrated radio spectrum of the Sausage 
relic seems to steepen at high frequenci es, perhaps more strongly than expected from ra¬ 
diative cooling alone fjStroe et al.ll2014b|) . So in the third model, we consider a power-law 
population with exponential cutoff as follows: 


fe,up (yP) 




(7) 


where 7e,c«t is the cutoff Lorentz factor. This may represent fossil electrons that have cooled 
down to ~ 7e,cut from the power-law distribution in Equation (|6]) . The integrated sp ectrum 
of the Sausage relic shows a strong curvature above 1.5 GHz fIStroe et al.ll2014bll . which 
corresponds to the characteristic frequency of synchrotron emission from electrons with 7 e ~ 
1.5 X 10^ when the magnetic held strength ranges 5 — 7 /iG (see Equation (IT^ below). So 
7e,cnt = 10^ and 2 x 10^ are chosen for SClpexl and SClpex2 models, respectively. 


Figure 3 shows the electron distributions in the cloud medium: the thermal distribution 
for the IGM gas with kT = 3.35 keV and the fossil population. The characteristics of 
the different models, SAl (hducial model), SAlp (power-law), and SClpexl (power-law 
with an exponential cutoff) will be described in Section 3.1. Note that there could be 
‘suprathermal’ distribution between the thermal and GR populations. But it does not need 
to be specihed here, since only /(p) around pinj controls the spectrum of re-accelerated 
electrons. 


3. RESULTS OF DSA SIMULATIONS 
3.1. Models 

We consider several models whose characteristics are summarized in Table 1. For the 
hducial model, SAl, the upstream density is assumed to decrease as p^p = po{r/r, 
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while the upstream temperature is taken to be kTi = 3.35 keV. Since we do not concern 
about the absolute radio flux level in this study, po needs not to be specifled. The preshock 
magnetic held strength is Bi = 2.5 pG and the immediate postshock magnetic held strength 
is B 2 {t) ~ 6.7 — 6.3 pG during 60 Myr, while the downstream held, i?dn(n), is given as in 
equation (j5]). The initial shock speed is Ug^i = 3.3 x 10^ km s“^, corresponding to the sonic 
Mach number of Mg^i = 3.5 at the onset of simulation. In all models with pup cx the shock 
slows down as Ug oc and at the electron acceleration age of 60 Myr, Ug ~ 2.9x10^ km s“^ 
with Mg ^ 3.1 and ainj ~ 0.73. The fossil seeds electrons are assumed to have a delta- 
function-like distribution around 7e,inj ~ 30. 

In SAlb model, the upstream magnetic held strength is weaker with Bi = 0.25 pG 
and B 2 {t) k, 0.67 — 0.63 pG during 60 Myr. Otherwise, it is the same as the fiducial model. 
So the character ‘b’ in the model name denotes ‘weaker magnetic held’, compared to the 
fiducial model. Gomparison of SAl and SAlb models will be discussed in Section 3.4. 

In SAlp model, the fossil electrons have a power-law population of /e,up oc while 
the rest of the parameters are the same as those of SAl model. The character ‘p’ in the 
model name denotes a ‘power-law’ fossil population. 


In SA2 model, both the preshock and postshock magnetic held strengths are constant, 
i.e., Bi = 2 pG and B 2 = 7 pG, otherwise it is the same as SAl model. 

In SA3 model, the upstream gas density decreases as pup = po{r so the shock 
decelerates more slowly, compared to SAl model. Gonsidering this, the initial shock speed 
is set to Ug^i = 3.0 x 10^ km s“^ with Mg^i = 3.2. At the acceleration age of 60 Myr, the 
shock speed decreases to Ug ~ 2.8 x 10^ km s“^ corresponding to Mg ^ 3.0. 


In SA4 model, the upstream density is constant, so the shock decelerates approximately 
as Ug oc more quickly, compared to SAl model. The upstream and downstream mag¬ 

netic held strengths are also constant as in SA2 model. Figure 2 compares the profiles of the 
flow speed, u{r,t), and the magnetic held strength, B{r,t), in SAl and SA4 models. Note 
that although the shocks in SAl and SA4 models are not very strong with the initial Mach 
number j 3.5, they decelerate approxi mately as Ug oc and n., oc respectively , 
as in the self-similar solutions of blast waves fjOstriker fc McKeelll988l:lRvu fc Vishniadll99ll ). 
The shock speed decreases by ~ 12 — 15% during the electron acceleration age of 60 Myr in 
the two models. 


In SBl model, the preshock temperature, Ti, is lower by a factor of 1.5^, and so the 
initial shock speed Ug^i = 3.3 x 10^ km s“^ corresponds to Mg^i ^ 5.3. The shock speed 
Ug ~ 2.9 X 10^ km s“^ at the age of 60 Myr corresponds to Mg ^ 4.6, so the injection 
spectral index ctinj 0.6. The ‘SB’ shock is different from the ‘SA’ shock in terms of only 
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the sonic Mach number. 


In SClpexl and SClpex2 models^ the fossil electrons have ^e,up P ‘exp^ (pie/'^e^cut^ 
with the cutoff at 7e,cut = 10^ and 'ye, cut = 2 x 10"^, respectively. The character ‘pex’ in the 
model name denotes a ‘power-law with an exponential cutoff’. A slower initial shock with 
Us,i ~ 2.2 X 10^ km s“^ and Ms,i ~ 2.4 is chosen, so at the acceleration age of 80 Myr the 
shock slows down to Mg ^2.1 with Oinj 1.1. The ‘SC’ shock differs from the ‘SA’ shock in 
terms of the shock speed and the sonic Mach number. The integrated spectral index at high 
frequencies would be steep with A^, ^ 1.6, while Ay ^ 0.7 at low frequencies due to the flat 
fossil electron population. They are intended to be toy models that could reproduce the in¬ 
tegrated spectral indices. Ay ~ 0.7 for v = 0.1 —0.2 GHz and Ay ~ 1. 6 for v = 2.3 —16 G Hz, 
compatible with the observed curved spectrum of the Sausage relic flStroe et al.ll2014b|l . 


3.2. Radio Spectra and Indices 


The local synchrotron emissivity, jy{r), is calculated, using the electron distribution 
function, fe{r,p,t), and the magnetic held prohle, B{r,t)- Then, the radio intensity or 
surface brightness, ly, is calculated by integrating jy along lines-of-sight (LoSs). 

rhmax 

Iy{R) = 2 / Jy{r)dh. (8) 

Jo 

R is the distance behind the projected shock edge in the plane of the sky, as dehned in the 
Introduction, and h is the path length along LOSs; r, R, and h are related as = (r^ — 
RY + The extension angle is 'ip = 10° (see Section 2.5). Note that the radio hux density, 
Sy, can be obtained by convolving R with a telescope beam as Sy{R) ~ Iy{R)7i9i92{l + z)~^ , 
if the brightness distribution is broad compared to the beam size of 9i92- 


The volume-integrated synchrotron spectrum, Jy = J jy{r)dV, is calculated by integrat¬ 
ing jy over the entire downstream region with the assumed geometric structure described in 
Section 2.5. The spectral indices of the local emissivity, jy{r), and the integrated spectrum, 
Jy, are defined as follows: 


a. 


(r) = 


A. = 


dinjyjr) 

dlnu 
din Jy 

dlnu 


(9) 

( 10 ) 


As noted in the Introduction, unresolved radio observations usually report the spectral 
index, ‘ex’ at various frequencies, which is equivalent to Ay here. 


In the postshock region, the cutoff of the electron spectrum, fe{r,'je), decreases down¬ 
stream from the shock due to radiative losses, so the volume-integrated electron spectrum 
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steepens from 7“*^ to 7e above the break Lorentz factor (see Figure 5). At the electron 
acceleration age tp„ p, the break Lorentz factor can be estimated from the condition t^ge = Aad 
( Kang et ah 2012 1: 


7e,br ~ 10' 


( ^age \ 

lOOMyr J \BI + 


( 11 ) 


Hereafter, the postshock magnetic held strength, il2, and -Brad are expressed in units of pG. 
Since the synchrotron emission from mono-energetic electrons with 7e peaks around 


/ 86-82 \ 2 


i/peak ~ 0.3 J 7 e ~ 0.63 GHz • 

the break frequency that corresponds to 7e,br becomes 


B, 




z/br ~ 0.63GHz 


/ Gge A 

VlOOMyry \Bl +Bl 


5 / VIOV ’ 


5 


( 12 ) 


(13) 


So the volume-integrated synchrotron spectrum, has a spectral break, or more precisely 
a gradual increase of the spectral index approximately from Oinj to Oinj -l- 0.5 around z/br. 


3.3. Width of Radio Shocks 


For electrons with 7e > 10“^, the radiative cooling time in Equation (j3]) becomes shorter 
than the acceleration age if tage ^ 100 Myr. Then, the width of the spatial distribution of 
those high-energy electrons downstream of the shock becomes 


A/cool(7e) ~ M2Gad(7e) ~ 100 kpc 


U2 \ 
10^ km s“^ / 



where U 2 is the downstream flow speed. 


(14) 


With the characteristic frequency Vpeak of electrons with 7e in Equation ([12]), the width 
of the synchrotron emitting region behind the shock at the observation frequency of z/obs = 
t'peak/(l + z) would be similar to A/cooi(7e): 


Al ^ 


W-Alcooihe) ~ IFTOOkpc- 


U2 


10^ km s 


-1 


52 + Bl 


rad 


fB2Y^" 

Vobs(l + z)' 

V 5 ) 

0.63GHz 


-1/2 


(15) 


Here, IF is a numerical factor of ~ 1.2 — 1.3. This factor takes account for the fact that 
the spatial distribution of synchrotron emission at z/peak is somewhat broader than that of 
electrons with the c orresponding 7e, because more abundant, lower energy electrons also 
make contributions flKang||2015al) . One can estimate two possible values of the postshock 
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magnetic field strength, 82 , from this relation, if Al^ can be determined from the observed 


prohle of surface brightness and U 2 is known. For example, Ivan Weeren et ahl fl2010l) inferred 
two possible values of the postshock magnetic held strength of the Sausage relic, i?iow ~ 
1.2 /iG and i?high ~ 5 /xG, by assuming that the FWHM of surface brightness, A/sb, is the 


same as ^ 


A/cooi(7e) ~ 55 kpc (i.e., W = 1). 


On the other hand, iKang et ahl (1201211 showed that the prohle of surface brightness 
depends strongly on 'ip in the case of the assumed geometrical structure. Figure 4 illustrates 
the results of a p/anar shock with diherent The shock has Ug = 2.7x 10^ km s“^ and Mg = 
4.5. The magnetic held strengths are assumed to be i?i = 2 /iG and B 2 = 7 /iG, upstream 
and downstream of the shock, respectively. The results shown are at the acceleration age of 
80 Myr. Note that the quantities, (/e(d, 7e), ju{d), and Iu{R), are multiplied by arbitrary 
powers of 7 e and Uohs, so they can be shown for diherent values of 7 e and i^obs with single 
linear scales in Figure 4. Here, the FWHM of (/e(7e) is, for instance, A/cooi(7e) ~ 28 kpc for 
7e = 8.47 X 10^ (red dotted line in the top panel), while the FWHM of ju is ^ 35 kpc 

for Uohs = 0.594 GHz (red dotted line in the middle panel). Note that the values of 7e and 
i^peak = ^'obs(l + z) are chosen from the sets of discrete simulation bins, so they satisfy only 
approximately the Equation (1121) . The bottom panel demonstrates that the FWHM of R, 
AZsb, strongly depends on -0. For instance, AZsb ~ 47 kpc for z/obs = 0.594 GHz, if -0 = 10° 
(red dotted line). This implies that due to the projection ehect, A/sb could be substantially 
diherent from A/cooi- So A/sb of the observed radio hux prohle may not be used to infer 
possible values of B 2 in radio relics, unless the geometrical structure of the shock is known 
and so the projection ehect can be modeled accordingly. 


Note that the quantity Y(B 2 ,z) = [(i?| + B^^)/5'^]~^ ■ {B 2 / 5 YR in Equation (IT^ also 
appears as in Equation ([13]). So the break frequency becomes identical for the two 
values of B 2 , By^w and i?high, that give the same value of Y. For example, SAl model with 
B 2 Y) ^ 6.7 —6.3 /iG and SAlb model with B 2 Y) k. 0.67—0.63 /xG (see Table 1) produce not 
only the similar AR^^^ but also the similar spectral break in Jy. But the corresponding values 
of 7e,br in Equation (1TT|) are diherent for the two models with Hiow and Hhigh- Moreover, the 
amplitude of the integrated spectrum would scale as Jy(i?high)/'^!y(-Biow) ~ (-Bhigh/Aow)^, h 
the electron spectrum ne(7e) for 7e < 7e,br is similar for the two models . We compare these 
two models in detail in the next section. 


3.4. Comparison of the Two Models with Hhigh and By^^ 

In Figure 5, we hrst compare the electron spectra and the synchrotron emission spectra 
in SAl and SAlb models. Here, the plotted quantities, //e, Gg, jy and Jy, are in arbitrary 
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units. The postshock magnetic held strength decreases from B 2 ~ 6.7 /iG to 6.2 fiG in 110 
Myr in SAl model, while B 2 ~ 0.67 /iG to 0.62 /iG in SAlb model. In these two models, 
the values of Y{B 2 ,z) are similar, so the spectral break in the integrated spectra should be 
similar as well. 

The left panels of Figure 5 show that the electron spectrum at the shock, ge{rs,p), steep¬ 
ens in time as the shock weakens. As expected, the volume-integrated electron spectrum, 
Ge{p) = p^Feip) = p‘^ f fe(r,p)dV, steepens by one power of p above the break momentum, 
Pe,br, which decreases in time due to radiative losses. In both models, however, the slopes, 
q = —d\nfe{rs,p)/d\np and Q = —dinFe{p)/dlnp, deviate from the simple steepening of 
Q = q to q + 1 above Pe,br, because of the time-dependence of the shock properties. 

The right panels of Figure 5 show the synchrotron spectrum at the shock, the 

volume-integrated synchrotron spectrum, J^, and their spectral indices, Q;inj(z^) and A^. Note 
that in both models, the transition of Ai, from ainj to ctinj -1-0.5 is broader than the transition 
of Q from g to g -|- 1. This is because the synchrotron emission at a given frequency comes 
from electrons with a somewhat broad range of 7e’s. As in the case of Q, A^ does not 
follow the simple steepening, but rather shows nonlinear features due to the evolving shock 
properties. This implies that the simple relation. Ay zz ctinj-1-0.5, which is commonly adopted 
in order to conhrm the DSA origin of synchrotron spectrum, should be applied with some 
caution in the case of evolving shocks. 

The highest momentum for ge(G,p), Peq oc Ug ■ [i?i/(i?g^ -|- 5^2)]^'^^) is higher in SAl 
model with stronger magnetic held, but the amplitude of ge{rs,p) near peq, say at pjmeC ~ 
10® — lO'^'®, is greater in SAlb model with weaker magnetic held. As a consequence, the 
ratio of the synchrotron emission of the two models is somewhat less than the ratio of the 
magnetic energy density, (i?high/-Biow)^ = 100. In SAlb model, for example, the amplitudes 
of jyir's) and Jy at 0.1 — 10 GHz are reduced by a factor of ~ 60, compared to the respective 
values in SAl model. Also the cutoh frequencies for both jy{rs) and Jy are lower in SAlb 
model, compared to those in SAl model. As pointed above, z/br is almost the same in the two 
models, although pe,br is diherent. This conhrms that we would get two possible solutions 
for B 2 , if we attempt to estimate the magnetic held strength from the value of z/br in fhe 
integrated spectrum. 


3.5. Surface Brightness Profile 

As noted in Section 3.3, the width of the downstream electron distribution behind the 
shock is determined by the advection length, A/adv ~ U 2 ■ tage, for low-energy electrons or by 
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the cooling length, A/cooi(7e) ~ ^2 ■ ^rad(7e), for high-energy electrons. As a resnlt, at high 
frequencies the width of the synchrotron emission region, varies with the downstream 

magnetic held strength for given U2 and i^obs as in Equation flT^ . In addition, the surface 
brightness prohle, Iu{R), and its FWHM, AZsb, also depend on the extension angle, ^|J, as 
demonstrated in Figure 4. 

In Figures 6 and 7, the spatial prohles of Iu{R) are shown for eight models listed in 
Table 1. Here, we choose the observation frequencies Uohs = 150, 600, 1400 MHz; the source 
frequencies are given as u = i^obs(l + z) with z = 0.192. While the downstream how speed 
ranges U 2 ~ 1, 000 — 900 km s“^ in most of the models, U 2 ~ 850 — 750 km s“^ in SClpexl 
model. So the results are shown at tage = 30, 60, and 110 Myr, except in SClpexl model 
for which the results are shown at tage = 30, 80, and 126 Myr. Note that the quantity, 
uR X X, is shown in order to be plotted with one linear scale, where X is the numerical 
scale factor specihed in each panel. In SBl model with a stronger shock, for example, the 
electron acceleration is more efficient by a factor of about 10, compared to other models, so 
the prohles shown are reduced by similar factors. 

For a hxed value of ip, the surface brightness prohle is determined by the distribution 
of ji/(r) behind the shock along the path length h, as given in Equation ([8]). The intensity 
increases gradually from the shock position (i? = 0) to the hrst inflection point, Rmipfl) = 
rs{t){l — cosip) 21 — 24 kpc with ip = 10°, mainly due to the increase of the path length 
along LoSs. Since the path length starts to decrease beyond if the emissivity, j^, 

is constant or decreases downstream of the shock, then R should decrease for R > 

In all the models considered here, however, at low frequencies, jy increases downstream of 
the shock, because the model spherical shock is faster and stronger at earlier times. As 
a result, R at 150 MHz is almost constant or decrease very slowly beyond or even 

increases downstream away from the shock in some cases {e.g., SA2 and SA4 models). So 
the downstream profile of the synchrotron radiation at low frequencies emitted by uncooled, 
low-energy electrons could reveal some information about the shock dynamics, providing 
that the downstream magnetic field strength is known. 

The second inflection in R{R) occurs roughly at i?inf,2 ~ W ■ AZadv for low frequencies, 
and at i?inf,2 ~ W ■ A/cooi(7e) for high frequencies where trad < fage, with W ^ 1.2 — 1.3. 
Here, A/adv ~ 100 kpc(M2/10^ km s“^)(tage/100 Myr), and AZcooi(7e) is given in Equation 
O- Figures 6 and 7 exhibit that at 30 Myr (tage < Aad), the second inflection appears 
at the same position (A/adv ~ 30 kpc) for the three frequencies shown. At later times, the 
position of the second inflection depends on tage at low frequencies, while it varies with B2 
and z/obs in addition to U 2 . Thus, only if tage > frad(7e), A/sb at high frequencies can be used 
to infer B2, providing that ip and U2 are known. 
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The width of the Sausage relic, dehne d as the FWHM of the surface brightness, was 
estimated to be AZsb rNj 55 kpc at 600 MHz (Ivan Weeren et allboioh . As shown in Figures 6 
and 7, in all the models except SA4, the calculated FWHM of Iu{R) at 600 MHz (red dotted 
lines) at the acceleration age of 60 — 110 Myr would be compatible with the observed value. 
In SA4 model, the intensity prohle at 600 MHz seems a bit too broad to £t the observed 
profile. 


3.6. Volume-Integrated and Postshock Spectra 


Figures 8 and 9 show the volume-integrated synchrotron spectrum, J^,, and its slope, 
A,^, for the same models at the same ages as in Figures 6 and 7. Here J^, is integrated over 
the curved, ribbon-like downstream volume with ip = 10°, as described in Section 3.2. In the 
figures, Jy is plotted in arbitrary units. The hlled circles represent the data points for the 
integrated flux of the Sausage relic, which were taken from Table 1 of IStroe et al.l (j2014b|) 
and re-scaled roughly to £t by eye the spectrum of SAl model at 60 Myr (the red d o tted line 
in the top-left panel of Figure 8). Observational errors given by IStroe et al.l (j2014b|) 
are about 10 %, so the error bars are in fact smaller than the size of the filled 
circles in the figures, except the one at 16 GHz with 25 %. 


At hrst it looks that in most of the models, the calculated reasonably hts the observed 
data including the point at 16 GHz. But careful inspections indicate that our models fail to 
reproduce the rather abrupt increase in the spectral curvature at ~ 1.5 GHz. The calculated 
Jy is either too steep at low frequencies < 1.5 GHz, or too flat at high frequencies > 1.5 
GHz. For instance, Jy of the fiducial model, SAl, is too steep at z/ < 1.5 GHz. On the 
other hand, SAlp model with a flatter fossil population and SBl model with a flatter 
injection spectrum (due to higher Ms) seem to have the downstream electron spectra too 
flat to explain the observed data for z/ > 1.5 GHz. 


As mentioned before, the transition of the integrated spectral index from Uinj to Uinj-1-0.5 
occurs gradually over a broad range of frequency, (0.1 — 10)z/br- If we estimate z/br as the 
frequency at which the gradient, dAy/du, has the largest value, it is ~ 4, 1, 0.3 GHz at 
Gge = 30, 60, and 110 Myr, respectively, for all the models except SClpexl model (for 
which Ay is shown at different epochs). Since z/br is determined mainly by the magnetic 
held strength and the acceleration age, it does not sensitively depend on other details of the 
models. In SClpexl model, the power-law portion (oc gives a hatter spectrum with 
Ay K, A at low frequencies, while the newly injected population at the shock with Mg ~ 2.1 
results in a steeper spectrum with Ay k. 1.6 at high f requen cies. Note that the spectral 
index estimated with the observed hux in IStroe et al.l fl2014b|) between 2.3 and 16 GHz is 
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Ay ^ 1.62, implying 2.1. 

In Figure 10, J^s from all the models considered here are compared at tage = 60 Myr 
(SA and SB shocks) or 80 Myr (SC shocks). The observed data points of the Sausage relic 
are also shown. In most of the models (SA models), the shocks have Mg ~ 3.0 — 3.1 at 
60 Myr, so the predicted Jy between 2.3 and 16 GHz is a bit flatter with Ay ^ 1.25 than 
the observed spectrum with Ay ^ 1.62 at the same frequency range. As mentioned above, 
SAlp and SBl models produce Jy, which is signihcantly flatter at high frequencies than 
the observed spectrum. The toy model, SClpexl, seems to produce the best fit to the 
observed spectrum, as noted before. 

In short, the volume-integrated synchrotron spectra calculated here steepen gradually 
over (0.1 — 10)z/br with the break frequency i^br ~ 1 GHz, if f^ge ~ 60 — 80 Myr. However, all 
the models considered here seem to have some difficulties fitting the sharp curvature around 
~ 1.5 GHz in the observed integrated spectrum. This implies that the shock dynamics 
and/or the downstream magnetic field profile could be different from what we consider here. 
Perhaps some additional physics that can introduce a feature in the electron energy spectrum, 
stronger than the ‘g -1-1’ steepening due to radiative cooling, might be necessary. 


In Figure 11, we present for the three models, SAl, SAlb, and SClpexl, the mean 


intensity spectrum in the downstream regions of [Ri, Ri -|- 5kpc], (R) = Iy{R)dR, 

where Ri = 5 kpc • (2z — 1) and i runs from 1 to 6. This is designed to be compared with th e 


‘spatially resolved’ spectrum behind the shock in radio observations (e.g.. lStroe et al.ll2013l) . 


The figure shows how the downstream spectrum cuts off at progressively lower frequencies 
due to radiative cooling as the observed position moves further away from the shock. 


4. SUMMARY 

We propose a model that may explain some characteristics of giant radio relics: the 
relative rareness, uniform surface brightness along the length of thin arc-like radio structure, 
and spectral curvature in the integrated radio spectrum over ~ (0.1 —10) GHz. In the model, 
a spherical shock encounters an elongated cloud of the IGM thermal gas that is permeated 
by enhanced magnetic fields and an additional population of fossil relativistic electrons. As 
a result of the shock passage, the fossil electrons are re-accelerated to radio-emitting energies 
(7e ^ 10^), resulting in a birth of a giant radio relic. 

In order to explore this scenario, we have performed time-dependent, DSA simulations 
of spherical shocks with the parameters relevant for the Sausage radio relic in cluster GIZA 
J2242.8-1-5301. In the hducial model, the shock decelerates from Ug ~ 3.3 x 10^ km s“^ 
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(Ms 3.5) to Us ~ 2.9 X 10^ km s“^ (Ms ~ 3.1) at the acceleration age of 60 Myr. The seed, 
fossil electrons with 7e,inj ~ 30 are assumed to be injected into the CR population, which 
is subsequently re-accelerated to higher energies. Such shocks are expected to produce the 
electron energy spectrum, fe{p) oc resulting in the synchrotron radiation spectrum with 

the injection index, ainj ~ 0.75, and the integrated index, A^, ^ 1.25, at high frequencies 
(> 1 GHz). We consider various models with a range of shock parameters, different upstream 
gas density prohles, different downstream magnetic held prohles, and three types of fossil 
electron populations, as summarized in Table 1. Adopting a ribbon-like curved shock surface 
and the associated downstream volume, which are con strained by the extension angle (or 


viewi ng depth) oi'ip = 10° as detailed in Section 2.5 fe.g lvan Weeren et al.ll2010l:iKang et ah 


20121) . the radio surface brightness prohle, A(R), and the volume-integrated spectrum, 


are calculated. 


The main results are summarized as follows. 


1) Two observables, the break frequency in the integrated synchrotron spectrum, i/br, 

and the width of the synchrotron emission region behind the shock, can have identical 

values for two values of postshock magnetic held strength (see Equations [13] and na). 

2) The observed width of the surface brightness projected onto the sky plane, AZsb, 
strongly depends on the assumed value of ip (see Figure 4). So A/sb may not be used to 
estimate the postshock magnetic held strength, unless the projection ehects can be modeled 
properly. 


3) The integrated synchrotron spectrum is expected to have a spectral curvature that 
runs over a broad range of frequency, typically for (0.1 — lOji^br- For a shock of Mg ~ 3 with 
the postshock magnetic held strength. Blow ~ 0.62 pG or Rhigh ~ 6.2 pG, the integrated 
spectral index increases gradually from Oinj ~ 0.75 to Oinj -l- 0.5 ~ 1.25 over 0.1 — 10 GHz, if 
the duration of the shock acceleration is ~ 60 Myr. 


4) Assuming that the upstream sound speed is Cg^i ~ 920 km s“^ (fcTi 3.35 keV) as 
inferred from X-ray observation, a shock of Mg ~ 3 and ~ 3 x 10^ km s~^ SA l 

model) can reasonably explain the observed width, A/sb ^ 55 kpc (Ivan Weeren et al.ll2010[) . 
and the curved integrated spectrum of the Sausage relic f Stroe et ah 2014bl) . SBl model 
with a shock of Mg ^ 4.5, however, produces the integrated spectrum that seems too hat to 
explain the observed spectrum above ~ 1 GHz. 


5) We also consider two toy models with power-law electron populations with exponen¬ 
tial cutohs at 7e 10^/e,up(p) OCp ^•2exp[-(7e/7e,c«t)^] (SClpexl and SClpex2 models). 
They may represent the electron populations that were produced earlier and then have cooled 
down to 7e ~ 10^. SClpexl model with a weaker shock [Mg k. 2.1) reproduces better the 
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characteristics of the observed integrated spectrum. But the steepening of the integrated 
spectrum due to radiative cooling alone may not explain the strong spectral curvature above 
1.5 GHz toward 16 GHz. 

6) This strong curvature at ~ 1.5 GHz may imply that the downstream electron energy 
spectrum is influenced by some additional physical processes other than radiative losses, 
because the integrated spectrum of radiatively cooled electrons steepens with the frequency 
only gradually. This conclusion is likely to remain unchanged even in the case where the 
observed spectrum consists of the synchrotron emission from multiple shocks with different 
Mach numbers, as long as the postshock electrons experience only simple radiative cooling. 

Other models that may explain the curved spectrum will be further explored and pre¬ 
sented elsewhere. 
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Table 1. Parameters for DSA Simulations 


Model 

Name 

Pup(r) 

Bi “ 

(/iG) 

B2{t) 

(/iG) 

kTi 

keV 

Us{t) 

(10^ km 

b 

Ms{t) ' 

/e,up (p) 

SAl 

_9 

oc r ^ 

2.5 

6.7-6.3 

3.35 

3.3 - 

2.9 

3.5- 

3.1 

— 

SAlb 

_9 

(xr 

0.25 

0.67-0.63 

3.35 

3.3- 

2.9 

3.5- 

3.1 

— 

SAlp 

_9 

xr 

2.5 

6.7-6.3 

3.35 

3.3 - 

2.9 

3.5- 

3.1 

oc 

SA2‘= 

_9 

xr 

2.0 

7 /iG 

3.35 

3.3 - 

2.9 

3.5- 

3.1 

— 

SA3 

-4 

xr 

2.5 

6.7-6.3 

3.35 

3.0 - 

2.8 

3.2- 

3.0 

— 

SAd'^ 

Po 

2.0 

7/iG 

3.35 

3.3 - 

2.8 

3.5- 

2.9 

— 

SBl 

_9 

xr 

2.5 

7.5 - 7.3 

1.49 

3.3 - 

2.9 

5.3- 

4.6 

— 

SClpexl'^' 

_9 

xr 

2.5 

5.5-5.1 

3.35 

2.2 - 

1.9 

2.4- 

2.1 

OC 

SClpex2^ 

_9 

xr 

2.5 

5.5-5.1 

3.35 

2.2 - 

1.9 

2.4- 

2.1 

OC p-^-"^ ■ ec2 


^Bi and B 2 are magnetic field strengths immediately upstream and downstream of the 
shock, respectively. See Equation (jl]). 

'’Since the shock is decelerating, the ranges of values are given for the acceleration age 
4ge = 0-60 Myr. 

’’The downstream (r < r^) magnetic field strength, B^n = B 2 in SA2 and SA4 models, 
while B^^{r,t) = B2{t) ■ [Pg{r,t) / Pg^2it)]^^'^ in the rest of the models. 

^£,1 = exp[-(7e/10^)'^]. 

^£,2 = exp[-(7e/2x 10^)2]. 
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Fig. 1.— Schematic diagram illustrating a scenario for the Sausage relic in a merging cluster, 
CIZA J2242.8+5301: a shock impinges on a fossil electron cloud with regular magnetic held 
in the cluster outskirt. Two gray disks represent the galaxy distribution, while the red 
ellipse shows the distribution of X-ray emitting gas. Arcs are meant to depict abundant 
shocks with different Mach numbers that are expected to form during a merger event. 
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Fig. 2.— Time evolution of the flow speed, u{r,t), and the magnetic held, B{r,t), of the 
spherical shock in SAl (top two panels) and SA4 (bottom two panels) models at the ac¬ 
celeration age, tage = 30, 60, 110 Myr (black solid, red dotted, and blue dashed lines, 
respectively). 
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Fig. 3.— Electron distributions in the cloud medium: the Maxwellian distribution of ther¬ 
mal electrons with kT = 3.35 keV (red solid line), fossil electrons with 7 e ~ 30 in SAl 
model (black solid), a power-law of in SAlp model (blue solid), and a power-law with 
an exponential cutoff in SClpexl model (green solid). See Table 1 for the different model 
parameters. The vertical dashed line demarcates the injection momentum, pinj. An unspec- 
ihed suprathermal distribution is shown in the red dotted line, but the electron distribution 
below pinj is not relevant for DSA. 
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0 20 40 60 80 100 


d, R (kpc) 

Fig. 4.— Spatial profiles of ge{d), ju{d), and Iu{R) for a steady planar shock with Ug = 
2.7 X 10^ km s“^ and Mg = 4.5 at the acceleration age of 80 Myr. Here, d is the downstream 
distance away from the shock, while R is the distant behind the shock projected in the sky 
plane. Top: 'y^'‘^^ge{d) is plotted for the electron Lorentz factor, 7 e = 5.98 x 10^ (black 
solid line), 8.47 x 10^ (red dotted), 1.29 x 10"^ (bine dashed), 1.59 x 10^ (green long dashed), 
4.24 X 10"^ (magenta dot-dashed), and 6.02 x 10^ (cyan dot-long dashed). Middle: I'^'^juid) 
is plotted for the observation freqnency, i^obs = 0.334 GHZ (black solid line), 0.594 GHz 
(red dotted), 1.33 GHz (bine dashed), 2.36 GHz (green long dashed), 16.7 GHz (magenta 
dot-dashed), and 29.8 GHz (cyan dot-long dashed). Bottom: u^-^R^R) is plotted for the 
extension angle ■0 = 5°, 10°, and 15° (solid, dotted, and dashed lines, respectively). The 
red (bine) cnrves are for Uohs = 0.594 GHz (1.33 GHz). The redshift of the host clnster is 
assnmed to be x = 0.192 here. 
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Fig. 5.— Comparison of SAl model with Bi = 2.5 pG (upper four panels) and SAlb 
model with Bi = 0.25 pG (lower four panels) at fage = 30, 60, 110 Myr (solid, dot¬ 
ted, and dashed lines, respectively). Left: Electron distribution function at the shock 
position, ge{rs,p) = p‘^fe{rs,p) (black lines), volume-integrated electron distribution func¬ 
tion, Ge{p) = p^Feip) = p"^ J fe{r,p)dV (red lines), and slopes of electron distribu¬ 
tion functions, q = —din fe{rs)/dInp (black lines) and Q = —dlnFg/dlnp (red lines). 
Right: Synchrotron spectrum at the shock position, i'jp{rs) (black lines), volume-integrated 
synchrotron spectrum, = v j jv{'f)dV (red lines), and synchrotron spectral indices, 
Oinj = —d\njy{rs)/d\ni> (black lines) and = —dlnJ^/dlnu (red lines). 
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Fig. 6.— Surface brightness profiles at i^obs = 150, 600, and 1400 MHz (black solid, red 
dotted, blue dashed lines) in SAl, SAlb, SAlp, and SClpex models (from top to bottom 
panels). See Table 1 for the model parameters. The results are shown at the acceleration 
age, fage, specified in each panel. For the extension angle, -0 = 10° is adopted. The quantity 
vlu X A is plotted with a scale factor X to present it in the same ordinate scale for all the 
models. 
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Fig. 7.— Surface brightness profiles at z/ = 150, 600, and 1400 MHz (black solid, red dotted, 
blue dashed lines) in SA2, SA3, SA4, and SBl models (from top to bottom panels). See 
Table 1 for the model parameters. The results are shown at tage = 30, 60, and 110 Myr (from 
left to right panels). For the extension angle, 'ip = 10° is adopted. The quantity x A is 
plotted with a scale factor X to present it in the same ordinate scale for all the models. 
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Fig. 8.— Volume-integrated synchrotron spectrum, Jj,, and its spectral index, at tage = 
30, 60, and 110 Myr (black solid, red dotted, and blue dashed lines, respectively) for SAl, 
SAlb, and SAlp models (top three panels), and at tage = 30, 80, and 126 Myr (black 
solid, red dotted, and blue dash ed lines, respectively ) for SClpexl model (bottom panel). 
Filled circles show the data from IStroe et al.l fl2014bf) . scaled to £t by eye the red dotted line 
(^age = 60 Myr) for SAl model. Observational errors are small, about 10%, except for the 
data at 16 GHz with 25% (shown in a vertical bar). 
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Fig. 9.— Volume-integrated synchrotron spectrum, Jy, and its spectral index, Ay, at tage = 
30, 60, and 110 Myr (black solid, red dotted, and blue dashed l ines, respect i vely) fo r SA2, 
SA3, SA4, and SBl models. Filled circles show the data from IStroe et al.l fj2014b|) . scaled 
as in Figure 8. 
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Fig. 10.— Volume-integrated synchrotron spectrum, J^, for all the models in Table 1. In 
the left panel, SAl (black solid line), SAlb (red dotted), SAlp (blue dashed), SClpexl 
(green long dashed), and SClpex2 (magenta dot-dashed) models are shown. In the right 
panel, SA2 (black solid line), SA3 (red dotted), SA4 (blue dashed), and SBl (green long 
dashed) models are shown. The results are shown at 60 Myr, except for SClpexl and 
SCl pex2 model s for wh ich the results are shown at 80 Myr. Filled circles show the data 
from IStroe et al.l fj2014bl ). scaled as in Figure 8. Observational errors (vertical bars) are 
small, about 10%, except for the data at 16 GHz with 25%. 
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Fig. 11.— Mean intensity spectrum, {1^), averaged over + 5 kpc] behind the shock, 

where Ri = 5 kpc ■ {2i — 1), i = 1, 2, 3, 4, 5, 6 (black solid, red dotted, blue dashed, green 
long dashed, magenta dot-dashed, cyan dot-long dashed lines, respectively) for SAl, SAlb 
and SClpexl models. The results are shown at the acceleration age, tage, specihed in each 
panel. 



















